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Abstract 

Explicit applications of factorization theorems for processes at hadron colliders near the hadronic 
endpoint have largely focused on simple final states with either no jets (e.g., Drell-Yan) or one 
inclusive jet (e.g., deep inelastic scattering and prompt photon production). Factorization for the 
former type of process gives rise to a soft function that depends on timelike momenta, whereas 
the soft function for the latter type depends on null momenta. We derive in soft-collinear effective 
theory a factorization theorem that allows for an arbitrary number of jets, where the jets are 
defined with respect to a jet algorithm, together with any number of non-strongly interacting 
particles. We find the soft function in general depends on the null components of the soft momenta 
inside the jets and on a timelike component of the soft momentum outside of the jets. This 
generalizes and interpolates between the soft functions for the cases of no jets and one inclusive 
jet. We verify consistency of our factorization theorem to 0(a s ) for any number of jets. While in 
this paper we demonstrate consistency only near the hadronic endpoint, we keep the kinematics 
general enough (in particular allowing for nonzero boost) to allow for an extension to partonic 
threshold resummation away from the hadronic endpoint. 
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1. INTRODUCTION 



Factorization of cross sections is the basis of every theoretical prediction at hadron col- 
liders. In its simplest form, factorization states the measured hadronic cross section a can 
be obtained by convolving a perturbatively calculable cross section a with nonperturbative 
parton distribution functions (PDFs) [H, 0], 

o = j®f®a. (1) 

The PDFs are universal, and can therefore be extracted from one process and used to make 
predictions in another. Moreover, a will in general depend on a hard scale Q (for example, 
the partonic center-of-mass energy yfs), while the PDFs depend on the scale at which they 
are measured, say Qq. The evolution of the PDFs between these two scales resums logarithms 
of Q/Qq. This basic paradigm illustrates the two main uses of factorization: separation of 
universal, nonperturbative contributions to a cross section from perturbatively calculable 
contributions, and resummation of logarithms of ratios of scales to which each contribution 
is sensitive. 

When a depends only on a single scale, Eq. fl2]) is the end of the story. The situation is 
more involved when a itself depends on multiple, widely disparate scales. For example, in 
many collider physics processes involving jets, a can depend on mass scales associated with 
the jets such as Mj, hard scales like y§, and seesaw scales like Mj/y§. In such cases large 
logarithms of ratios of these scales can spoil the convergence of the fixed order perturbative 
expansion of a. One must further factorize a in order to resum these large logarithms 
and, perhaps, to separate out any other nonperturbative physics that is not captured in the 
PDFs 0-0]. The general structure of resummation, at the level of next-to- leading logarithms, 
has been previously investigated in Q. 

In this paper we will focus on so-called threshold logarithms. When a process approaches 
its kinematical threshold, there is limited phase space available for radiation. This gives 
rise to an incomplete cancellation between real and virtual diagrams, resulting in large 
logarithmic terms. This is common in situations where the invariant mass of the final 
state is near the maximum available energy, which limits the amount of energy that can 
go into excess radiation. Examples of this type of resummation can be found for Drell- 
Yan, deep-inelastic scattering (DIS), B meson decay, and top production 0, [s|-24]. It has 



been suggested [25|, [26j that a similar effect occurs at hadron colliders away from hadronic 



endpoint due to the steepness of parton luminosities, and this effect was explored more 



quantitatively in [27|. In this paper, we will concentrate on hadronic threshold and assume 
that the invariant mass of the final state is near the maximum allowed by the collider; 
however, we plan on exploring resummation away from the hadronic endpoint in future 



work [28]. For this reason, we will derive a factorization theorem that can be applied away 
from hadronic threshold without loss of information. 

An extremely useful tool to prove factorization is effective field theory. In the case of 



jet physics, soft-collinear effective theory (SCET) [2l|, |29M31| is the relevant effective field 
theory that can be used to derive factorization in many hard scattering processes 32) . The 
SCET Lagrangian is constructed by integrating out all modes of QCD except for soft modes 
and collinear modes with respect to some fixed number of directions ?v Matching QCD 
onto SCET gives rise to a hard function that contains the physics of the hard scales in the 
problem, and matrix elements of the remaining soft and collinear fields give rise to soft and 
jet functions, respectively. 
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The first applications of SCET involved cases with particularly simple jet definitions, such 



as in Drell-Yan [2 



or 



33 where there are no jets, hemisphere jets in event shapes [34T4441 
completely inclusive jets as in B — > X s ^ 2jJ, DIS |22|-|24| . and prompt photon production 
[45l | . Factorization of jets defined with more generic algorithms was considered in (ifl and 
two-jet rates defined with jet algorithms were computed using SCET in (35[ l47l. ]48|. A 
study of various different jet algorithms and the dependence on the jet parameters in the 
framework of SCET was discussed in [48]. More recently, an NLL analysis of jet shapes in 



multijet events using modern jet algorithms in e + e collisions was performed in |49l . 150 



The goal of this paper is to derive a factorization formula for an arbitrary number of jets 
in the presence of any number of non-strongly interacting particles in the threshold limit. 
We allow for a nonzero total rapidity and calculate the ingredients of this formula to allow 
resummation of threshold logarithms at NLL accuracy. This is conceptually distinct from 
the case of a single final state jet, which can be measured indirectly by simply demanding 
that a non-strongly interacting particle is produced with nonzero pp- When there is more 
than one final state jet, jet algorithms must be used to identify jets, and so the technology 
of incorporating jet algorithms into a factorization formula, developed in 46[ and applied 



to e + e collisions in 49, 50 



must be employed. The consistency of this factorization 
(that is, the fact that the cross section is independent of the factorization scale /1) is only 
demonstrated here in the true hadronic endpoint. However, we plan on investigating the 
consistency of this factorization away from hadronic endpoint using the steepness of parton 



luminosities in 28 



In this paper, we will assume that Glauber modes do not contribute at leading order in 
the power counting. This cancellation has only been formally proven at the level of cross 



sections for simple processes, e.g. |5l|. While Glauber modes could potentially spoil our 



factorization theorem, we assume that for sufficiently inclusive observables the argument in 



5l| generalizes. We also assume that PDFs can always be factorized from the partonic cross 
section, i.e. that 

da da 

= f f . (2) 

dO J J dO K J 

is always valid. This has also not been proven to be generically true, but is phenomenologi- 
cally valid for a wide range of observables. 

Factorization formulas for the case of a single, inclusive jet have soft functions that depend 
on the null component of the total soft momentum in the direction of the jet. When there 
are no jets (e.g., Drell-Yan), the soft function depends on the timelike component of the 
total soft momentum. In extending threshold resummation to more than one jet using jet 
algorithms, we find a soft function that depends on the timelike component of the total soft 
momentum outside of the jets and on the null component of the soft momentum within 
each of the jets. Thus, our result reduces to the previously considered cases of zero and 
one inclusive jet when our jet algorithm either includes none or all of the final state soft 
momentum, respectively. 

The organization of this paper is as follows. In Sec. [21 we define precisely what we mean 
by threshold production of N jets and discuss the corresponding kinematics. In Sec. [21 we 
discuss how the definition of threshold affects the logarithmic structure of the result and 
how, in particular, our definition should not introduce so-called non-global logarithms. In 
Sec. HI we briefly discuss different classes of jet algorithms used at hadron colliders. We then 
derive our factorization theorem in Sec. [51 beginning for notational simplicity with the case 
of a single (quark or gluon) final state jet, then extending these results to the case of N jets. 
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We derive the anomalous dimensions for the objects that appear in our iV-jet factorization 
formula in Sec. |6]and use these in Sec. [7] to show that our factorization theorem is formally 
consistent, at least in the hadronic endpoint region. Finally, we present our conclusions in 
Sec. El 



2. KINEMATICS OF THRESHOLD RESUMMATION 

To explain our approach to threshold resummation, how it includes both the cases of 
Drell-Yan and direct gauge boson production as limiting cases, and how it is extendable 
to arbitrary iV-jet production, we first discuss the kinematics. By demanding that the 
final state contains N jets each with fixed transverse momentum p^ and pseudorapidity 77, 
together with some number of non-strongly interacting particles with total 4-momentum q, 
we are requiring that there is a minimum partonic center-of-mass energy 

N 



where p l j is the momentum of the ith jet. This momentum is defined in terms of the and 
i] j of the jet as 

Pj = (p J cosh r?j,pj, p J sinh 77 j) . (4) 

Of course, the actual partonic center-of-mass energy s typically exceeds this minimum value, 
and in general can be as large as the available machine center-of-mass energy s. Therefore, 
the dimensionless variable z, defined as 

•^min / r \ 

can range from 

t < z < 1 , with r = . (6) 

s 

Going to hadronic threshold (r — >■ 1) forces z — > 1, such that the only emissions kine- 
matically allowed are collinear radiation off the hard partons that form jets, as well as soft 
radiation. Radiation collinear to one of the jets with momentum scaling as E cm (l, A 2 , A) 
(in the light-cone coordinates of the jet) and soft radiation scaling as E cra (X 2 , A 2 , A 2 ), each 
contribute an equal amount to s, where A ~ a/1 — -2 is a small, dimensionless parameter. In 
this limit of restricted radiation, partonic momentum conservation can be written as 

Pl = q+ k s +j2pi, (7) 

i 

where pj is the total initial-state (partonic) momentum, k s is the total soft momentum, and 
p % c is the momentum carried by collinear fields in the direction of jet i. The total momentum 
can be separated into two components: the first is the minimum momentum needed to 
create iV jets of fixed pr and 77 together with the non-strongly interacting particles of total 
momentum q, while the second brings the invariant mass of the final state above its minimum 
value ,§ m i n . To do this, we note that an arbitrary four-vector p can be written as the sum 
of a massless four-vector which characterizes the transverse momentum and pseudorapidity 
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of p and a purely timelike four- vector with a magnitude equal to the +- component of p in 
light-cone coordinates about n = (l,p/|p|) (i.e., p + = p° — |p|), which characterizes the 
off-shellness of p. That is, for any four- vector p, we can write 

p» = p> 1 J + p+v* , (8) 

where v M = (1,0) and pj is given in Eq. with p"j and rjj the transverse momentum and 
pseudorapidity of p, respectively. We want to apply this relation to the total 4-momentum 
in each of the jets. To do this, we note that the jet algorithm will group some of the soft 
momentum k s into parts that belong to jet i, fcj, and a part that is not included in any of 
the jets, fc out , 

N 

k s = S ^k i + fc out . (9) 

i 

Using this together with the relation Eq. (jSJ) and the fact that 

(p« + h) + = pf + fc° - |p« + k<| = pl + + kf + 0(A 4 ) , (10) 

where on the right-hand side, kf is plus with respect to pj and p* + is plus with respect to 
p c , we can write momentum conservation Eq. ([7]) at leading order in A as 

N N 

rf = 5" + Ct + J2p" + v "[J2(pt + k t)] • (ii) 

i i 

Here, we have also used that out-of-jet collinear radiation is power suppressed jj50j]. 
Given these definitions, we can write 

1 - z = - §P i ■ (fcout + v [ J2(pt + H)] ) + 0(\ A ) 

N 

• (vi ■ fcout + P? + k t)) + > (I 2 ) 



s 

where s = pj. We see that since pi is timelike, 1 — z depends on the timelike component of 
the soft momentum outside of the jets and on the null component of the momentum within 
the jets. 

So far, we have discussed the kinematics in the hadronic endpoint defined as r — > 1, which 
is the main focus of this paper. However, z can be forced close to 1 not only in this hadronic 
endpoint, but also in the limit of steeply falling parton luminosities. In this case, final states 
with small values of s are preferred, giving again z — > 1. Our analysis is independent of the 
precise mechanism which guarantees that 1 — z can be regarded as a small quantity, and can 
therefore be used away from the true hadronic endpoint. 

We are now in a position to discuss how our parametrization of 1 — z reduces to the 
standard variable in the case of Drell-Yan and cases when there is one inclusive jet, such 
as B — » X s 7, DIS, and direct gauge boson production. As we will see in Sec. El Eq. (TT2j) 
implies the soft function in general depends on the timelike component of fc out and on the 
null components (with respect to the corresponding jet directions) of the soft momenta in 
each of the jets fcj. In Drell-Yan, there are no jets in the final state and so the entire soft 
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momentum k s is just k out . This is why the soft function in Drell-Yan depends only on the 
timelike component of the total soft momentum. For a single inclusive jet (i.e., defined 
with a jet algorithm that includes all of the hadronic momentum), all the soft momentum is 
included in the jet, such that k s = k%. This explains why the soft function in this latter type 
of process only depends on the null component of the total soft momentum. In Ref. p], on 
the other hand, threshold resummation for dijet production was considered and it was found 
that the soft function only depended on the timelike component of momentum outside of the 
jets. This apparent discrepancy is due to the fact that the limit of small jet size R — > was 
taken and the contribution of in-jet soft particles vanishes in this limit. 1 From the discussion 
above, the soft function will have dependence on the null component of in-jet momentum 
for jets of finite size. 



3. NON-GLOBAL LOGARITHMS AND THE DEFINITION OF THRESHOLD 

In Ref. jsj, two definitions of z were defined for the case of two final state jets, which we 
denote here and Zb, 



Z„ 



s 



z b = ; , (13) 

s 

where pi j2 are the total 4-momenta of the jets . Note that, unlike pj (cf. Eq. (j3j)), p 12 can 
not be defined by the jet direction and energy alone, so both definitions of z are indirectly 
sensitive to the jet mass. To lowest order, 1 — z a ^ can be written as 

?k° 

1 _ ^o»t + 0(A 4) 



M 



j j 

'-- 2 -k + ^m7 +0{x>) ' (14) 

where on the right-hand side, Mjj can be set to {pi + P2) 2 for both z a> b to order A 4 . 

It is well-known that jet observables which have an energy scale of radiation inside of a 
let that is widely disparate from the scale outside the jet give rise to non-global logarithms 



52j. This is due to the fact that in this case, real emission corrections to radiation inside 



the jet can be vetoed when one of the daughter particles escapes the jet, and this leads to 
an incomplete cancellation of real and virtual corrections, an effect which is stronger for 
radiation near the jet boundary (the so-called "buffer region"). As has been pointed out in 
the literature (see, e.g., Ref. |53j), z a introduces non-global logarithms, whereas Zb does not. 
This is clear because the limit 1 - 4 < 1 only restricts radiation outside of the jet and, 
for a jet size R ~ 1, the radiation within the jet is not restricted in the hadronic endpoint 
r — > 1. On the other hand, the scaling of in-jet and out-of-jet radiation is correlated with 



1 Note that double counting is avoided in [6( by removing collinear modes from the soft function ( "eikonal 
subtractions"), which, for R ~ A, removes any dependence on the soft momentum inside of jets, at leading 
order in the power counting. 
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zj, (in particular, in both cases soft radiation has momentum components that scale as 
Mjj(1 — Zf,)) such that no non-global logarithms should arise. 

Notice that for a two-jet final state in the hadronic endpoint (r — > 1), the definition of z 
given in Eq. 03) reduces to 

1 - Z —^Mj] + 2 ^ + 0{X) > (15) 

where we used Eq. (jT^J) and that, to C(A 2 ), p/ = (V§,0) = (Mjj, 0) in the r — >• 1 limit. 
Thus, our definition of z restricts radiation both inside and outside of the jet, similarly to 
Zb in this limit and should correspondingly not introduce non-global logarithms for R ~ 1. 
(This is also clear directly from Eq. ( fl2|) which is valid away from the hadronic endpoint, 
provided another mechanism enforces 1 - 2 C 1 such as the steepness of PDFs.) 



4. JET ALGORITHMS AT HADRON COLLIDERS 

Perturbative calculations require a precise definition of the phase space boundaries im- 
posed by the jet algorithms. Any jet algorithms needs to be infrared safe; otherwise, it 
leads to infrared divergent results when calculated in perturbation theory. There are two 
general types of jet algorithms: cone algorithms and cluster algorithms. Cone algorithms 
decide on which particles belong to a given jet based on cones of fixed size R, while cluster 
algorithms group particles together into jets based on a relative measure of their distance. 
These jet algorithms act on the entire set of particles in the final state to decide how many 
jets are contained in a given event and which particles belong to which jet. Almost all jet 
algorithms depend on a jet size R, and a distance AlZij that measures the distance between 
two particles in r\ — cf> space 

Ally = yj (Ar^jj) 2 + (A^ij) 2 . (16) 

As already discussed in the previous section, the relevant degrees of freedom in jet pro- 
duction close to z — 1 are collinear and soft particles. To perform perturbative calculations 
in this region we therefore need a restriction on these degrees of freedom to decide whether 
they belong to a given jet or not. Collinear particles in a given direction all belong to the 
same jet. This is in contrast with soft particles, which can either belong to a jet or not. 
Note that the treatment of jet algorithms in SCET is only correct to leading order in the 
power counting parameter A. Therefore, we assume that all jets have energy much in excess 
of their mass, and that all jets are widely separated. 

Standard cone algorithms, such as SISCone, are quite simple. The restrictions they 
impose on each particle to belong to a given jet are independent of other particles in the 
event, and only depend on the angular distance from the jet direction. The restriction for 
both soft and collinear particles to be in a jet j with direction n is therefore 

©s R of tJ = ©fen,, = n e ( A7 ^« < ® ■ 

i 

For the purposes of this paper, we only need results at relative order a s , and therefore only 
have to consider one extra particle in the final state. The restrictions therefore simplify, and 
for the extra particle i we can write 

©X- = ©foil,, = ©( A ^,n < R) ■ ( 18 ) 
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Cluster algorithms iterate a process of calculating a distance measure d^ for all pairs of 
particles and a distance d{ for each jet and, if a di is smallest, removing particle i, or, if a dij is 
smallest, merging jets i and j. In general , d^ is defined as dy = min{di, dj}ATZij/D, where 
D is a fixed parameter that characterizes the size of a jet and the precise definition of di 
depends on the choice of algorithm. 2 This makes the action of the jet algorithm considerably 
more complicated. As explained above, all collinear particles in a given direction have to 
end up in the same jet, which allows us to write a generic restriction for the action of a 
cluster algorithm on a set of collinear particles in a given direction as 

AT-l 

©foil. = n e ( A7 ^ fe < R) • (19) 

k=0 

Here, A denotes the total number of collinear particles in the direction of the jet rij. ATZ^ k 
denotes the AlZij between the pair of collinear particles in the set of N — k remaining 
particles with the smallest d^. For soft particles, such a generic formula is not possible 
(at least analytically), since different soft particles can end up in different jets, and the 
restriction on a given particle depends on all other soft particles in the event. At relative 
order a s , however, the restrictions 0f for cluster algorithms simplify and are given by 

6f oftJ = 0(A7^„ < R) (20) 

for soft particles and by 

ef 0llJ = e(AK k>l < r) (21) 

for collinear particles, where k and I label the new particles after the collinear splitting. For 
more information about jet algorithms in SCET, see [50| . 



5. iV-JET FACTORIZATION THEOREM 

In this section, we present the factorization theorem for the cross section to produce 
A jets, defined with respect to a jet algorithm, differential in the 3-momentum (pp and 
pseudorapidity rj) of each jet and of the non-strongly interacting particles. To keep the 
notation simple, we begin in Sec. 15 Al by discussing the case of a single jet produced via 
the channel qq — > g. We then discuss the differences between this derivation and the one 
needed for the channel qg — > q. It will be clear from these derivations that, aside from the 
promotion of the hard and soft functions to matrices which arise from mixing of operators in 
color space, there is nothing conceptually or technically new for arbitrary iV-jet production 
in our approach. This allows us to generalize our results to the A-jet factorization formula 
in Sec. 15 Bl 

In writing down a factorization theorem, we first assume that we can match QCD onto 
operators in SCET containing N + 2 distinct collinear fields. This is valid when a (direct or 
indirect) measurement constrains the final state to be A-jet like. In our case, this is ensured 
by the fact that we take the variable 1 — z to be small, together with the assumption that 
the jets are well separated from each other and from the beams (with the latter requirement 



2 An example for such a distance measure is di = pl^ for the fcx algorithm. 
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ensuring that the probability of initial-state collinear radiation to produce a jet is power 
suppressed relative to the probability of the jet arising from the hard interaction). Our 
derivation is agnostic as to the cause of 1 — z <C 1, and in [28] we explore in greater detail in 
what regimes the steepness of parton luminosities allow the factorization theorem we derive 
here to be applied away from the hadronic endpoint. We will assume in this section that 
the reader has some familiarity with SCET. For details, we refer the reader to the original 
SCET literature HH, MM. ■ 



A. Case of a Single Jet 

1. qq^g 

Working to leading order in the electroweak coupling constant, we first write the full 
theory matrix element mediating the partonic interaction as 

(qXlOlP^) = J2 M^T a A b (Xl^A^P,) . (22) 

% 

Here, \q) represents the non-strongly interacting final state of total momentum q, \Pi) and 
IP2) are the incoming hadrons with the corresponding momentum, and \X) represents the 
hadronic final state. The index i labels Dirac structures. This equation defines the M°^ M . 
Note that we have used the fact that there is only one color singlet in the decomposition of 
3® 3® 8. 

In terms of M i; the matching of QCD onto the fields of SCET takes the form 

=Y. M f"Y. c ^y { * i+ ^ > ( 23 ) 

3 {f>} 

where pi is the label momentum carried by the field. At tree level, we have 

C*;({P}) = V (24) 
The matching condition in momentum space takes the form 

Mf"Q a ^(k) = Mf" J d A xe- ik - x Q a ^(x) 
j ^i=i J / 

x (2n) 4 5 A (p 1 +p 2 -p 3 -k), (25) 

where we turned the sums over labels and integrals over residual momenta into integrals 
over the full d 4 pi and used the shorthand notation 
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Using this matching condition, we can write the cross section differential in the jet pseu- 
dorapidity r]j and transverse momentum as 



da 



rest. 



d 2 pjdtanh rjjd§ q 2E 2 a 



Y,\(qx\o\p x p 2 



I spin avg 



6{ VJ - r)(X))5*(p§ - p T (X)) 



(27) 



cm x 
4r4/ 



x (2n) A 5\P 1 + P 2 -q- Px 

_ rest. 

—EE M f" 



2EL 



x 



-__ rpArpA 



k=l 



spin 



x 5( VJ - v(X))5 2 (pj - p^(X)) (2 7 r) 4 5 4 (p 1 + p 2 - p 3 - q) . 

Here, we defined d$ g as the phase space measure of the m non-strongly interacting final-state 
particles, 



n 



2Eu 



(28) 



where we have omitted symmetry factors coming from identical particles in $ g . The re- 
striction on the sum over final states X ("rest.") is that they include exactly one jet as 
defined by the jet algorithm and the delta functions that fix rj(X) and p T (X) act on the 
part of X identified to be the jet, which we assume to be sufficiently separated from the 
beams such that contributions from collinear initial-state radiation are power suppressed. 
We also define M = 7°M t 7°. To arrive at this equation, we used the hadronic momentum- 
conserving delta function in the first line to shift the operator to the point x, applied 
the matching condition Eq. (1251) . and then integrated over x which resulted in the partonic 
momentum-conserving delta function on the last line. 

We can simplify Eq. f[2T|) using the following observations. First, we can use the BPS 
field redefinition [3l( to decouple soft and collinear modes to 0(X 2 ), 



Xn{x) -> Y n (x)Xn{x) 
Xn(x) -> X n (x)Y 7 l(x) 

B n (x) -> y n (x)B n (x) = Y r \{x)B n {x)Y n {x) , 
where Y n is a soft Wilson line, for which we adopt the conventions 3 



P exp 



ds n ■ A s (x + sn) 



for an incoming particle (where P denotes path ordering) and 



x 



Pexp 



\g s I ds n ■ A s (x + s n) 
Jo 



(29a) 
(29b) 
(29c) 



(30a) 



(30b) 



3 For a discussion of the various conventions for incoming and outgoing Wilson lines and how they are 
related see, for example, Ref. 54 1. There is also a nice discussion of how soft Wilson lines arise in the 
path integral formulation of SCET in Appendix C of Ref. 27 1. 
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for an outgoing particle. In momentum space, the field redefinition in Eq. (129)) induces the 
shift p c — )■ p c + k s , where k s is the total soft momentum carried by the Wilson lines. 

Second, we can implement the restriction on the final state X that there is one jet at the 
operator level to all orders by using the jet algorithm operator 9f , defined in Sec. HI This 
allows us to complete the sum over X and factorize the collinear and soft matrix elements 
from one another. 

The final state collinear matrix element can be written as 

(0\Bl A (-p' 3 ) Q^AVJ - V)S 2 (P T J - P T ) B* (- ft )|0) = -Ej(2n)W(p' 3 - p,)6 AB g^J°M , 

(31) 



which defines the gluon jet function J^{p). Integrating Eq. f l3Tj) over the full p% and p' 3 , 
which contain integrals over residual momentum and sums over the labels u and n 3 of the 
B± field, fixes the labels to be n% = (l,pj/|pj| coshr/jjtanhr/j) and u = 2p^ cosh rjj. The 
label "3" on 0^ and on indicates these are defined with respect to the direction n 3 of 
the jet. Note that this jet function depends on the choice of the jet algorithm. 

For an algorithm that is inclusive over collinear initial-state radiation, the initial-state 



collinear matrix elements give rise to PDFs from the relations |32|, |46 



iW{p\r a/ {p')x:{p)\p)^^. 

fpfp'(P\X a a'(-p)X a a(- P )\P)s 



F A it"*" 



/spin avg. 



dxf q (x) 



E rm S nn >{-) I dxf q (x) , 



(32) 



with px 2 set to Zi^-Ecm-^ = |^i,2 n \,2 (where rii=n and n 2 = n) and p[ 2 = p lj2 . 

The color structure in the matrix elements Eqs. (I3T1) and (I3"2"|) leads to a trace over the 
color structure of the Wilson lines in the soft function. We can write the soft function in 
terms of the variables of interest k ont and k 3 (the out-of-jet and in-jet momenta) as 



^k s S(k s , {rii}) 



d L(0 T 



yt yt T A V y 

1 n 2 1 n 3 J - 1 n 3 1 n 1 



CaCf 

"^out d hS(k ont ,h,{ni}) 



yt yt t a V y 

1 n 1 1 n 3 ± 1 n 3 1 n 2 



(33) 



where {rii} = {n q , rig, n g } and T (T) denotes (anti-) time ordering. S(k ont , k%, {rii}) is then 
defined as 



S(k out ,k 3 ,{ni}) 



CaCf 



-<0|T 



yt yt T A V y 

1 n 2 1 n 3 ± 1 n 3 1 n\ 



S 4 (k 3 - h) 



X 5 4 (/Cout - ^out)T 



yt yt t a V y 

i m «3 1 n 3 1 n 2 



|o> 



(34) 



and it is understood that the replacement k s — > k out + k 3 should be made wherever k s 
appears. The operators k out and k 3 are defined as 



^out = ( 1 ~~ @Sft,3 ) 



(35) 
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Finally, using that the spin- and color-averaged square of the Born matrix element Mb, 



\Mb\ 2 } can be written in terms of Mi as 



M. 



1 



B 



2 = s% £, |M 

spin, color 



B\ 



Tr[T A T M 



J2M^MT(pir(p^g 



3 



4C A ^ 



(36) 



we see that the Dirac structure of the matrix elements in Eqs. fl3Tl) and fl32l) naturally gives 
rise to the Born cross section. 

To simplify the notation, we define a hard function H which includes all perturbative 
corrections contained in the matching coefficients Cy as 



H({n h u,i}) = 



Op 
4CV 



3 M f [ — 



\M f 



(37) 



where are the labels on the three collinear fields. H by definition is 1 to leading order in 
a s . Putting this together, we arrive at the expression 



da 



E, 



d 2 pjdtanh77jd$ g 2£ c 2 n 



/ ^^ |Mb|2 / i% J / d% (38) 



x H({m, uj i })f q (x 1 )fq(x 2 )J^(ps)S q q-, g (k out , k 3 , {rii}) 



n 



x (27r) 4 5 4 (x 1 E cm - + x 2 E, 



n 



q - Kut - Pa- h) 



The final step is to simplify the momentum-conserving delta function. We use Eqs. (jSJ) 
and (fTTj) to write it, up to power corrections in A, as 



n 



n 



S (xiE cm - + x 2 E cm - - q- k out -p 3 - k 3 ) 



2 ' 2 



2 Pl ■ k out + 2p%p+ + k^ 



x 5 2 (pj +qT)5(F-tanh- 1 



p J sinh Tjj + q z 
p T j cosh Tjj + g 



(39) 



where pj is the total (partonic) initial-state momentum, p® = (x% + x 2 )E cm /2 and p\ 
{x\ — x 2 )E cm /2. Here, we used that 



5[f(x,y)]6[g(x,y)] 



S(x - x )S(y - y ) 



d(f,g) 
d{x,y) 



(40) 



where f(x ,y ) = g{x ,y ) = and made the change of variables 



x 2 



—e 



Y 



—e 



(41) 
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In switching from to z and the total rapidity Y, we will also need that 

' dxi dx2 



Xi x 2 



l dz 

z 



dr. 



(42) 



Now, since there is no dependence on the soft momenta other than on the components 
Pi ■ k out and k% and the only unconstrained component of p% is the plus-component p£, we 
can integrate over the other variables to obtain our final expression 



da 



7T 



d 2 pjdtanh r]jd§ q 



E 4 T 

cm 



dz 



2 T 



dY \M B \ 2 H q ^gx({n i ,Ui})f q (x 1 )fq(x2) 



x JdptJ 9 M) |dCt/d4^ 9 (Ct^ 3 + )5 2 (pJ + q T ) 
.! / p T j sinh^j + q z 



X 8 \Y — tanh 



p j cosh Tjj + q 



x S 



l-z-^= (cosh Y (pi + kfi + Ct) 



(43) 



The soft function in Eq. 



is defined as 
1 



CaCf 



T 



V"f Vt T^V V 
"2 «3 - , n3" I «l 



5(fc 3 + -fc 3 + )MCt- 



x T 

<S(Ct)<W) + 0(a.) . 



Vt vt T^V V 

n 1 1 nz ± 1 n 3 1 n 2 



(44) 



Note that, despite our notation, fc° ut = • fc ut/|p/| is simply the projection of the soft 
momentum outside of all jets onto a timelike vector. This function is most easily computed in 
the partonic center-of-mass frame where pi is in fact purely timelike. However, in general pi 
will have nonzero spatial components, in which case k® ut is not simply the energy component 
of the k out four- vector. The jet function in Eq. (143]) is defined as 

JZ(p + )9T6 AB = ~^f d'xe^(0\B&(x)€^B»£(0)\0) = 5(p + )g?5 AB + 0(a s ) , (45) 



where, again, the label u is set to w = 2Ej = 2p^coshr]j. 



2. qg-^q 

The majority of the above discussion goes through in much the same way for the channel 
qg — > q. The main differences are that the matrix element of final state fields gives rise to a 
quark jet function, defined by 

/ i% (o|x£&) - v)s 2 (p t j - p t )x c M\o) = 5 C ' C (^) J*fa) , (46) 
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which, after integrating over all but the p + component, becomes 

J«(p + )S ab = ^J d 4 xe^-(O|| X ^Ja ; )0f olli 3^,jO)|O) = 5(p+)8 ab + 0{a s ) . (47) 
For the initial-state gluon, we obtain the gluon PDF via the relation 

J fpfp>(P\Bt A (p')B*(p)\P) = j^Sab jf 1 ^f g (x) . (48) 

Finally, the soft function is defined as in Eq. ( 133|) but with the appropriate modification in 
the definition of the Wilson lines for incoming and outgoing fields given in Eq. ( )30l) . The 
final result after these differences are taken into account is of the form Eq. (T4"51) but with 
the substitutions f q —> f g , J 9 — > J q , and S qq ^ g — > S qg ^ q . 



B. Extension to N Jets 

The above results clearly generalize. The only nontrivial complication is due to mixing in 
color space. To avoid cumbersome notation, we will explain what generalizes and state the 
result rather than write down the iV-jet derivation. Explicitly, for every quark, antiquark 
and gluon in the final (initial) state, the all-orders jet function (PDF) has precisely the Dirac 
and color structure to contract with the non-QCD matrix element M^ '^ " to give the Born 
matrix element at tree level. The PDFs and jet functions have the same definitions as in 
the single jet case. However, due to the fact that operators with different color structures in 
general mix, the hard and soft functions in these formulas should be interpreted as matrices 
in color space. 

Since the main difference in the generalization to N jets is in the soft function, we will 
discuss it in more detail. It takes the general form 

1 k N 

S(kl v {ni,kt}) = ^<0|O^(Ct - ?Ll*p*) JJ W _ kt)O 8 \0) , (49) 

where M is a normalization factor such that the soft function is unity (times delta functions 
in its arguments) at tree level and the operator 0$ is the product of Wilson lines that arise 
from the BPS field redefinitions Eq. (f29|) . appropriately traced. As for the one jet case, the 
arguments hi in the soft function run over all final state jets and the rii include all directions, 
both initial and final. The operators A; out and k{ that appear in Eq. (f49|) are defined as 

K = @soft,i ^ 
N 

&t= (i-E^ft,)^- (50) 

i 

The result of going through the same steps as for the single jet case leads to the iV-jet 
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factorization formula, 

da 

nf d 2 pj.dtanh r^d"^ 



E 4 T 

era 

N 



dz I dY\M B \ 2 H({n i ,u i })f 1 (x 1 )f2(x 2 ) 
dpi 



N 



2-(2vr) : 



MpDU hK dkl t s(kl t ,{ ni ,k+}) 



x (2tt) 4 5 



1 - z 



N 

(coshYj2(pt + kt)+k 



o 

out 



x 5 \ Y — tanh 



/ v^W _|_ • i . \ \ N 



(51) 



In Eq. ( 15 ip . both H and S are matrices in color space, while J and / are proportional to 
the identity in color space and there is an implicit trace. Loop corrections will in general 
mix color structures, which means that beyond tree level H and S will contain off diagonal 
elements. The details of resummation in the presence of color mixing is discussed in [7|. 

6. ANOMALOUS DIMENSIONS 

Now that we have shown that a generic iV-jet cross section factorizes in the limit of 
1 — z <C 1, we are left to calculate each ingredient of the factorization theorem. In this 
paper, we focus on the consistency of the factorization theorem to 0(a s ), and we therefore 
need the one-loop anomalous dimensions for the hard, jet and soft functions. Note that the 
results presented here are enough to resum threshold logarithms at NLL 4 . 



A. Hard, Jet, and Parton Distribution Functions 

Both the Born-level matrix element and the hard function, which can be found by cal- 
culating the virtual corrections to the Born-level matrix element in the MS scheme, are 
process dependent, so they can not be calculated generically. However, the hard anomalous 
dimension, defined as 



dH({ni,Ui};tJ,) 
dln/x 



(52) 



only depends on the directions m, label momenta Wj, and color charges T, of the collinear 
particles and is given by (|55j, 0, Hsj, @) 



j H ({ni, Ui}; fi)= ^2 (- r cus P T 4 2 In ^ - 7^ j - 2 r cusp ^ TV Tj In 

ie{partons} i / 



n j • rij 



(53) 



4 There are several different ways to define precisely what is meant by NLL. In this paper, we will use the 
convention of 27 1. 
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The cusp anomalous dimension at two loops is given by 



cusp 



a s at 
— H - 

71 ATX 2 



67 



71" 
~3 
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The 7, to one-loop are given by 



7, 



3a, 



C f 



and 



1 9 



9 



a\ 



2tt y 2tt 

for quarks and gluons, respectively, where (3q is defined as 



A) 



/So 



11C, 



2iV 



(54) 



(55) 



(56) 



The sum on is a sum over all distinct pairs of partons i and j for z ^ 7. 
The quark and gluon jet functions have been calculated previously, e.g.[60l- 

HI, 

and were 

first calculated with a jet algorithm in 50] . Their anomalous dimensions, defined by 



dJi(Po/x) 
dln/i 



(57) 



are given by 



7Ji (p+ ; a*) = (2 r cusp t 2 in — + 7 <) - 2 r cusp Tf i C 



1 -1 4 

p 4 



(58) 



The expressions for r cuS p and 7^ are the same as for the hard function. Note that the 



algorithm in [50j used a polar angle for the measure and not Eq. ( II 6p . However, since 
the anomalous dimension in that case did not depend on the algorithm parameter R, the 
result must be independent of which measure is chosen since the precise definition of the jet 
boundaries is not associated with any singularities. It does however affect the finite parts 
of the jet function, which become important starting at NNLL accuracy. 

It is well-known that the parton distribution functions are not perturbatively calculable; 
in practice, they can be expressed as universal matrix elements, which are then extracted 
from experiment. However, the evolution of the PDFs with fi can be computed, 



d/i(xi;/x) a s f 1 dz /a* 



(59) 



where the repeated index j is summed over and are the Altarelli-Parisi splitting functions. 
Near hadronic threshold, the splitting functions simplify and can be written as 



-Pij(x) 



2 r t 2 

" 1 cusp 1 j 



1-x) 



+ 7; 5(1 - X) 



(60) 



B. Soft Function 



In general, the soft function depends on the null component of the soft momentum inside 
each jet as well as the timelike component fc° ut (defined in Sec. [5] as k® nt = pi ■ k out /\pj\). At 
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order a s , the soft function can be written as a sum of functions that depend only on one 
momentum variable, with trivial dependence on the others and is given by 

s(ki t ,{n t ,kt})=s out (ki t ) n 5(kt)+5(k° j e s ^ k t) n ( 6i ) 

iS{jets} «e{jets} ie{jets} 

where the sum over i G {jets} is over all i corresponding to outgoing jets and does not 
include the incoming partons (and we remind the reader that the dependence here on ki 
is only over final state jets but the run over all initial and final partons). The timelike 
component of the soft function, S out , receives contributions from soft gluons that are not 
inside any of the outgoing jets. We can find this by calculating the contribution of soft 
gluons going anywhere and then subtracting the contribution from gluons that enter one of 
the jets. This can be written in the hadronic center-of-mass frame as 



5ou t (Ct) = -E T ^ T ^ 2 ^v e / 



d k rii-rij 



(2-n) d {m-k){n r k) 
x2vr^ 2 )5(^-Ct)[l- E © 



R 

soft,fc 



A;S{jets} 



(62) 



where 9f is the restriction that the gluon is in jet k, defined by a jet algorithm of size R as 
in Sec. SJ This is most easily calculated in the partonic center-of-mass frame. Denoting the 
directions and energies of the collinear partons in this frame as hi and tDj respectively, we 
have that 

go B *(0 = -E T «- T i^/^ f . 2v5(k>)S(k°-kl t )\l- e ©sV 

(63) 

where we have used the fact that, for an r\ — cf> algorithm, the jet algorithm restrictions are 
frame invariant. 

The null components of the soft function, S- m (kf), are defined in the hadronic center-of- 
mass frame as 

S in (K)--J2 T ^ 2 9 2 ^ I f ftf y h ^ (k 2 )5(n k .k-kt)Q^ k . (64) 

J (27r) d (ni-k){nj-k) 

In the partonic center-of-mass frame, this can be written as 
SUH) = -E T.T, 2 9 V I ^^-^2^)^ - e 5(a . (65) 

The calculation of S{ n and 5 0U t in the partonic center-of-mass frame can be related to 
the calculation of the soft function in jHo[. While [50| uses a polar angle measure, we will 
show that, even though Si n and S out separately depend on the parameter R, the anomalous 
dimension is R independent. This implies that all algorithms with the same singularity 
structure as the polar angle algorithm used in [50[ have the same anomalous dimension for 
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this observable. Specifically, the anomalous dimension calculated in the partonic center-of- 
mass frame should be the same for both a polar angle measure and an rj — measure. Since 
the r) — 4> measure is boost invariant, the R independence of the anomalous dimension must 
be true in all frames 5 . The relations of S in and S out to the soft function in |50j are given by 



Sout(^out) — 2 



(id) 



dA 



oincl 



fee {jets} 



(66) 



J A=Ct 



and 



(67) 



where Slj, S£j, and 5'™ eas are defined and computed in [50jj. In calculating Sfj eas , we use 
the definitions Tg = k^/tik and Sr = S (tq — k + /uk), where r a and 8r are originally defined 



in 



50 



The anomalous dimension of this soft function is defined as 



dS(k° ut ,{n h kt};y) 
dln/x 



"■out 



xS(kZ t ,{ni,k'l};n) 



(68) 



Using the results of |50j, together with Eqs. fl66l) and fl67|) . the result for the anomalous 
dimension can be written as 



7s(C t ,KA + }^)= E ^M;v) W + 7s oat (C t ;^) II 

iG{jets} j6{jets} i€{jets} 



(69) 



with 



7Si (/c+;/i) = 2r cusp T^ 



2 A*w< 



7s out (Ct! /*) = Tcusp E T r Ti 2 In ^ 5(Ct) " 4 r cusp (T? + T 



,2a 1 / M 

■2 



(70) 
(71) 



7. CONSISTENCY OF FACTORIZATION TO 0{a s ) 

Consistency is a nontrivial check of our factorization theorem. The factorized cross 
section should be independent of the factorization scale fi in the threshold limit and thus 
renormalization group invariant. Starting from the generic iV-jet cross section, Eq. fl5Tl) . 



5 We have verified this explicitly for small R by making the replacement R — > R/coshi], which relates the 
polar angle and rj — <p measures in the small R limit. 
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ignoring multiplicative factors that do not affect the derivative, and using the shorthand 
notation 

da d da 



d In /i d In /i & 2 Pj- drjj, d^ 



(72) 



we have that 
da 



dln/i 



oc / dz 



±ln - 

2 m T 



dY# ({ri;, Wi}; /i)/i(xi; fi)f 2 (x 2 ; n) 
II JdptJi(pt;fi) J] Jdkt Jdk out S(k out ,{n u kt}^) 



*S{jets} 



ie{jcts} ' 



dlnF({ri i ,a; i }; i w) dln/i(xi;/x) , dln/ 2 (^2;^) 



dln/i 



dln/i 



dln/x 



+ dlnJj(pt ; ^) ^ dlnff(fc° ut ,{^,fc+};/i) 



x 5 



f - z 



dln/i dln/t 
-^(coshF J2 (Pt + kt) + k 



iG{jets} 

2 



o 

out 



iG{jets} 



(73) 



There are several simplifications we can make to check the independence of /i. First, /i 
only enters perturbative expressions, and whether or not the cross section depends on /i is 
independent of nonperturbative physics. This allows us to use the perturbative definition 
of the parton distribution functions. Second, given that the /i dependence of each of the 
factorization ingredients starts at order a s , we can use the tree level expressions for the 
hard, jet and soft functions, as well as for the PDFs, 



/<(x;/i) = 5(1 - x) 
H{{n h Ui}]n) = 1 

J%(pt;n) = $(pt) 

S(k out ,{n l ,k+};ri = 5(kl t ) H *(*+). 

iS{jets} 

Using this and working to lowest order in a s , we can simplify Eq. (1731) to get 
da 



(74) 
(75) 
(76) 

(77) 



n ^Pn(r) + ^P 22 (T) + - 
d m /i 7r 7r 2 cosh Y 



E 

iG{jets} 



2 cosh Y 



(1 - r);/i 



+ 



2 cosh Y 



E 

i€{jets} 



Is, 



2 cosh Y 
1 - r);/i 



(l-r);/i + 1h {ii)5{1-t) 



(78) 



After plugging in Eqs. ( 1531 . ( 1581) . and ( )69|) . rescaling the plus functions using the identity 

1 



ax 



In a . 1/1 

S(x) + - - 

a a \ x 



(79) 
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and combining the various terms, we find 



da 

d OC 1 cusp 
m \i 



i , 2 i j 2 



YTt-Tj 21n^ +T?ln^+T 2 2 ln^ + V T? 21n- 



«6{jcts} 



After making the simplification 



rii ■ rij 



*,3 



a;,- u;. 



i€{partons} 



(80) 



(81) 



where we have used that pi-pj = PvPj and J2i = 0, Eq. ( J80l) gives 

da 



dln/i 



0. 



(82) 



This result confirms that our factorization theorem is consistent at hadronic threshold and 
justifies using the renormalization group to resum logarithms of 1 — r. 



8. CONCLUSIONS AND OUTLOOK 

We have derived a factorization theorem for the production of N jets, together with 
any number of non-strongly interacting particles, such as electroweak gauge bosons. This 
factorization theorem allows us to write the physical cross section in terms of a convolution 
of parton distribution functions, a hard function, and jet functions for each observed jet and 
a soft function describing among other things the color recombination between the initial 
and final state partons. Both the jet and the soft functions depend on the precise form of 
the jet algorithm chosen. 

The main new ingredient in this factorization theorem is a soft function that depends on 
a timelike component of the soft momentum outside of the observed jets, and the lightlike 
component of the soft momentum in a given jet. This function is directly related to the 
soft function first proposed and calculated for the case of jet production in e + e~ collisions 
in 49|, |50j. This soft function allows us to interpolate between the soft function arising for 
final states without observed jets (which depends only on a timelike component of the soft 
momentum) and the soft function for completely inclusive jet production (which depends 
only on the lightlike component of the soft momentum). 

We have derived the UV divergent parts of all ingredients of the factorization theorem 
to 0(a s ). These were then used to show that the combination of all ingredients of the 
factorization theorem is independent of the arbitrary factorization scale /x, and therefore 
the derived results satisfy the nontrivial requirement of consistency. While consistency was 
shown in this work only in the true hadronic endpoint, we have kept the kinematics general 
enough (in particular allowing for a nonzero overall boost) to allow for a generalization of 
our results to the case where the steepness of the parton luminosities force events to be close 
to the partonic threshold. This result, which is by far more interesting phenomenologically, 
will be the subject of future work (28| . 

Our results can be used to explicitly resum threshold logarithms to NLL accuracy (in 
the log-counting convention of [27]]) for any process in hadron collisions with any number of 
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jets and non-strongly interacting particles in the final state. The technology of going from 
the anomalous dimensions we present here to explicit resummed distributions is well-known. 
Beyond one jet, in addition to the standard resummation methods, we need the matrices 



Tj-Tj, but these have been computed for many processes (see, e.g., [6l. l55l. l63H66j ) . including 



all 2 — > 2 and 2 — > 3 partonic channels. The only other ingredient needed to obtain a NLL 
distribution is the Born matrix element. 

In addition, if our results are extended to include two-loop results of the anomalous 
dimensions together with the full algorithm- dependent one-loop finite parts, NNLL results 
can be obtained for all processes for which the virtual NLO corrections are known. Together 
with recent advances in calculations of NLO cross sections (e.g., W + W~j [67, 68|| and Wjjj 
(69|). this would have a significant impact on the precision frontier of predictions at the 
LHC. 
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